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I. INTRODUCTION 



A future molecular nanotechnology has to rely on techniques and systems that enable 
us to selectively interact with single molecules and prepare them in defined states. Only 
then molecules can act as well-defined building blocks for nanotechnological devices, e.g. for 
highly integrated data storagd^ or in photomechanical machined^HS! \ class of molecules 
that fulfills these requirements in solution are so-called molecular switches, i.e. molecules 
that can be reversibly switched between two or more stable states via external stimulP. In 
the quest to exploit this intriguing functionality for nanotechnological devices, research has 
recently concentrated on the properties of such molecules when they are localized at solid 
surfaced^. Herein, adsorption at metal surfaces represents a particularly appealing sub- 
topic, with couplings to the underlying substrate that go beyond the mere weak physisorption 
limit. Notwithstanding the constant danger of ultrafast quenching of any (photo-)excitation 
due to a too strong coupling, the very intricacies of this interaction with the delocalized 
metal electron manifold is hoped to give rise to novel isomerization mechanisms and thereby 
to a switching behavior that is not attainable in gas-phase or solution. 

A corresponding switching has for instance been achieved for azobenzene derivatives at 
coinage metal^^'^and the geometric structure of the (meta)stable surface mounted molecu- 
lar states is to some extent unravelecP^^. The observed non-trivial excitation wavelength, 
surface orient at iorP^H^ and functionalization!^^^ dependence of the switching efficiency in- 
deed suggests that the isomerization process proceeds either on molecular potential energy 
surfaces (PESs) that are strongly modified as compared to their established counterparts in 
solution, or even involves electronic excitations in the underlying metaP^. In this situ- 
ation independent insight as provided by material-specific first-principles theory would be 
highly desirable to further elucidate the mechanistic details of the switching function. Un- 
fortunately, the very experimental evidence already indicates the tremendous challenge that 
metal-surface mounted switches pose to such modeling: On the one hand, the theory ob- 
viously needs to accurately describe both molecular ground and involved excited electronic 
states. Particularly the latter is commonly the realm of numerically highly-demanding cor- 
related wave-function based approaches tractable only for very limited system sizes. On the 
other hand, the non-trivial infiuence of the substrate dictates its explicit treatment, which 
in order to properly describe the metal band structure needs to rely on extended supercell 
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geometries. Together with the sheer lateral extension of flat-lying adsorbed molecules like 
azobenzene this gives rise to system sizes that are already at the cutting-edge of what can 
be tackled with approximate ground-state techniques like density-functional theory (DFT) 
with present-day semi-local exchange-correlation (xc) functiona y^^ * ^^ * ^^ * ^^ * ^°l This calls for 
numerically highly efficient approaches to describe the excited states, which in fact should 
only impose CPU-costs comparable to those of a semi-local DFT ground-state calculation. 

It is self-evident that such approaches will be approximate in nature, and thus need to 
be carefully scrutinized to assess what can and what cannot be addressed reliably. A careful 
scrutiny requires accurate references as benchmark though. With only limited and indirect 
experimental information on excited PES topology available this primarily concerns higher- 
level theory. Notwithstanding, corresponding techniques might have their own limitations. 
As little conducive as uncritically applying approximate theories is then to readily dismiss 
them either for what they are or because of discrepancies with incorrect reference data. With 
this scope the objective of the present work is specifically to revisit the reliability of the nu- 
merically undemanding DFT-based A Self Consistent Field (ASCF) approach to excited 
stateP^'^in the context of the prototypical molecular switch azobenzene (Ab, H5C6-N=N- 
CgHs). Apart from perfectly meeting the computational efficiency requirement this type of 
constrained DFTf^ZEH technique is particularly appealing as it is readily extended to apphca- 
tions at extended surfaces, even in case of appreciable adsorbate-surface hybridizatiorP^'^. 
Unfortunately, a preceding study by Tiago et alW^ on gas-phase azobenzene reported quali- 
tatively different ASCF PES topologies as compared to Time-Dependent DFT (TD-DFT)13Z1 
reference data, questioning the usefulness of this approach for an envisioned application to 
the isomerization mechanism of azobenzene (and its derivatives) at metal surfaces. 

Revisiting the problem, we employ approximate Coupled Cluster Singles and Doubles 
(RI-CC2)E3I3i! as additional reference technique. Its accuracy in describing the lowest lying 
singlet excitations relevant for the isomerization of gas-phase azobenzene was already demon- 
strated by Fliegl et al^. For reasons of computational feasibility the present calculations are 
still exclusively performed for the free molecule in the gas-phase. Nevertheless, the discus- 
sion and assessment will also be geared towards the application of ASCF to surface mounted 
azobenzene. As such, we e.g. concentrate on semi-local DFT with a gradient-corrected xc 
functional (GGA)ESlas ground-state basis for the approach, as such functionals are presently 
the unbeaten workhorse for metal surface studies. From the mapping of different ground- 
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and excited-state PES cuts discussed in the context of azobenzene photoisomerization we 
arrive at the central conclusion that ASCF, TD-DFT and RI-CC2 yield overall very similar 
PES topologies. The discrepancies reported before by Tiago et al. were caused by an undis- 
covered state crossing in their TD-DFT calculations. Limitations due to the approximate 
xc functional in ASCF and TD-DFT primarily show up as global excited-state PES offsets, 
leading to a severe underestimation of absolute vertical excitation energies at e.g. the level 
of semi-local functionals. 

II. COMPUTATIONAL METHOD 

All ground-state calculations at the spin-polarized Kohn-Sham (KS) DFT level were 
performed with the all-electron full-potential DFT code FHI-aim^. Centrally targeted are 
results as obtained with the GGA functional due to Perdew, Burke and Ernzerhof (PBE)'3Slto 
describe electronic exchange and correlation. In order to assess the effect of different xc treat- 
ments additional calculations were performed with the local-density approximation (LDA) 
functional in the parameterization by Perdew and Wan^^ and with the hybrid functional 
g2LYpE2l40| FHI-aims employs basis sets consisting of atom-centered numerical orbitals, 
and we specifically used the "tier2" basis set with the internal default tight settings for the 
numerical integrations. From test calculations at the "tierS" basis set level we conclude that 
calculated relative energies (and the ensuing ASCF excitation energies) are converged to 
within 2meV at these numerical settings. The code also offers a standard implementation 
of the ASCF approachl^^'^ to obtain approximate excited states. The basis of this method 
are enforced modifications of the population of individual KS states, i.e. the KS equations 
are solved self consistently under the constraint of a given excited state occupation of the 
KS states, thereby accounting for orbital relaxation. To make full excited state PES scans 
and geometry optimizations tractable within ASCF we modified this discrete constraint to 
a Gaussian smeared constraint that affects the population of all KS states within a defined 
small energy window of 0.01-0.02 eV width. At generally insignificant changes of the total 
energy, this allows us to also readily converge systems with degenerate KS states that other- 
wise lead to significant problems in the SCF procedure. The actual geometry optimizations 
in both ground and excited states were then performed in combination with a locally modi- 
fied version of the Atomic Simulation Environment (ASE)^^. Forces were hereby relaxed to 
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below a threshold value of 10 me 

In the most straightforward ASCF realization, singlet excitations are mimicked by en- 
forced population changes within one spin channel, while triplet excitations are modeled 
through appropriate occupation changes in both spin channels. In terms of the azoben- 
zene frontier orbitals {vide infra), the singlet SI state would thus result from the enforced 
occupation of the KS lowest unoccupied molecular orbital (LUMO) and the enforced de- 
population of the KS highest occupied molecular orbital (HOMO) within one spin channel, 
while the triplet Tl state would result from the enforced occupation of the KS LUMO in 
one spin channel and the enforced depopulation of the KS HOMO within the other spin 
channel. However, earlier detailed workP^'^^HSl j^as demonstrated that the resulting single 
determinants of KS orbitals yield a particularly inaccurate description precisely of low-lying 
singlet states as are of central interest here. This owes to the fact that present-day local and 
semi-local xc functionals only evaluate the electron density that is symmetry-broken with 
respect to the true multi-determinantal singlet states. A possible remedy to this that we will 
employ throughout this work is the so-called "sum method" (SM) of Ziegler et al^, in which 
the multiplet corrected energy of the singlet state is calculated from the single-determinant 
singlet and triplet energies as = 2Es — -Et- A second possibility is to calculate a singlet 
state simply using non-spinpolarized DFT calculations, in which case the magnetization 
density is zero everywhere in spacd^. In such restricted DFT calculations the occupations 
of the involved doubly-degenerate KS states are then simply varied by ±1. While this ap- 
proach lacks a proper formal justification, it is particularly appealing in the context of the 
envisioned calculations for surface-mounted azobenzene, as there non-spinpolarized DFT 
would represent a significant saving in computational time. 

The TD-DFlSSand RI-CC2 calculation^™! ^sed to assess the performance of the ASCF 
approach were done with TURBOMOLE V6.2^^. TD-DFT calculations were hereby done 
for the same xc functionals as in the ASCF well as for the long-range corrected 

CAM-B3LYP functional^'', which is presently not available in FHI-aims nor in TURBO- 
MOLE. The CAM-B3LYP calculations were therefore performed with the GAMESS cod(P. 
In addition to these excited state calculations TURBOMOLE V6.2 was employed to obtain 
the constrained PBE geometries for the PESs in Section IIIB[ For these and the TD-DFT 



calculations we used a gaussian basis set of triple zeta quality (def2-TZVP) with polariza- 
tion functions from the Ahlrichs series of basis functionJ^, and the resolution of identity 




FIG. 1. Schematic overview of two possible photoisomerization mechanisms for azobenzene: 
rotation around the central CNNC dihedral angle to (upper panel) and inversion around one of the 
NNC angles a (lower panel). 

(RI) approximation!^. We estimate that the relative energies and excitation energies are 
converged within 10-20 meV with respect to a quadruple zeta basis set. For the TD-DFT 
calculations the maximum value of the eucledian norm of the residual vector for the tran- 
sition density matrices was set to 1-10~^. For the RI-CC2 calculations the basis set was 
def2-TZVPP, which was previously shown to yield highly accurate excitation energies for 
the azobenzene systempSl. With respect to the higher def2-QZVPP basis set we estimate the 
uncertainty in the calculated relative and excitation energies to be about 60 meV. 

III. RESULTS AND DISCUSSION 

A. Isomerization of gas-phase azobenzene: Current state-of-the-art 

In order to unequivocally set the target for the approximate ASCF approach we first 
briefly recapitulate the present state of understanding concerning the switching behavior of 
azobenzene and its derivatives, as it has emerged from a plethora of previous experimental 
and theoretical studies. This archetypical molecular switch is characterized by two stable 
isomers, namely cis (Z) and trans (E) azobenzene, whereby E-Ab is more stable by 0.6 eV 
and both states are separated by a sizable ground-state barrier of about 1.6 eV with respect 
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to the E-Ab state^. Both isomers are interconvertable via photoexcitatiorP^. Azobenzene 
can Z— T-E or E— t-Z isomerize following excitation in the UV-visible regime to either the so- 
called Sl(n —7- TT*) or the S2(7r — )■ tt*) state^^. The corresponding quantum yields for E— t-Z 
(Z— i-E) isomerization in n-hexane are 0.40 or 0.12 (0.53 or 0.25) following excitation to SI 
or S2, respectivei y^^*^^ l Recently it has been shown that excitation to higher lying states 
predominantly leads to dissociatiorP'. 

The detailed electronic and nuclear dynamics behind molecular Ab switching, espe- 
cially photoinduced switching, have undergone a vast amount of experimentaP^*^ and 
theoretica F^ * ^^ * ^^* ^ investigation over the past decades and are still to some extent controversiaP^EMMlIlIl^ 
Here, time resolved femtosecond absorption and fluorescence measurements pose particu- 
larly efficient tools to investigate the underlying electron dynamics from the experimental 
side. Braun and coworker^ measured transient absorption spectra and assigned several 
time constants to different parts of the isomerization process. For the SI (S2) isomerization 
of E-Ab the authors find: ri =0.34(0.42) ps, r2=3.0(2.9) ps, r3=12(12) ps. In the case of S2 
excitation it was necessary to assign an additional smaller time constant of 0.13 ps. From the 
similar time constants they concluded that both processes happen over fast nuclear motion 
from the Franck-Condon (FC) structure to the conical intersection (CI) (ri), relaxation to 
the ground state minimum energy structure {T2) and further vibrational cooling through 
the solvent (t^). The authors interpret the additional process following S2 excitation as 
fast initial population transfer from S2 to SI. In this view, both excitation channels thus 
follow the SI dynamics. For the isomerization starting from Z-Ab they reach the same 
conclusion even though the time constants show more subtle differences between SI and 
S2 (ri=0. 17(0.2) ps, r2=2(l.l) ps, r3=10(14) ps). These experimental values have been 
confirmed by several other groups even though the corresponding authors came to different 
interpretations concerning the actual isomerization mechanismP^. From their calculations 
of the involved PESs, Cattaneo and Persico^^ as well as Ishikawa and coworker^ also pro- 
posed that excitation to S2 is immediately followed by transitions that could involve several 
states and finally reach the SI state, from whereon the actual isomerization then follows SI 
dynamics. Recently this has been further supported by high level multireference calculations 
of Conti, Garavelli and Orlandi^^ and explicit non-adiabatic dynamics simulations by Yuan 
and coworker^. 

For the nuclear SI dynamics several possible pathways have been discussed in the litera- 
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FIG. 2. Two-dimensional relaxed ASCF(GGA-PBE) PES scans of rotation around the dihedral 
CNNC angle (jJ and inversion around one of the two CNN angles a, cf. Fig. [TJ Shown are the 
ground state SO (left), the first excited state SI (center) and the second excited state S2 (right). 
Energies relative to the zero reference E-Ab ground state energy are given in eV. 

ture. Figure [T] illustrates the two most frequently studied mechanisms, namely an isomeriza- 
tion around the central CNNC dihedral angle ("rotational pathway") and an isomerization 
around one of the two CNN angles ("inversion pathway"). The general understanding of the 
prevalence of these mechanisms has undergone various transitions. The initial belief was that 
excitation to SI mainly follows inversion whereas excitation to S2 should follow rotational 
isomerizat ionl^^'^, this way rationalizing the different quantum yields. However, most of the 
recent experimental studied^^HlS]^ .^g^ theoretical studies that either investigated the 
excited PESj^^ ° * ^^ * ^^ l or performed explicit non-adiabatic dynamics simulation d^^ * ^^ * ^^ * ^^ * ^^ * ^^ ! 
agree on the dominance of the rotational isomerization following excitation in either SI or 
S2 for azobenzene in gas-phase and solvent. 

B. Ground and excited state PES topology 

As established in preceding detailed quantum chemical work the centrally targeted low- 
lying excitations, SI and S2, have largely singly excited character at least at the FC struc- 
tures, and can be viewed as n ^ vr* (HOMO to LUMO) and vr ^ vr* (HOMO-1 to LUMO) 
transitions, respectivelj^^^^. In SM-corrected ASCF we accordingly model the two singlet 
excitations by modifying the populations of HOMO and LUMO (SI) and HOMO-1 and 
LUMO (S2), respectively. Figure [2] shows correspondingly obtained two-dimensional PES 
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scans along the dihedral CNNC angle u and along one of the two CNN angles a, cf. Fig. [T} 
Each point of the PES corresponds to a molecular geometry, in which the values for these 
two angles were constrained to the specific value, while all other degrees of freedom of the 
azobenzene molecule are those as resulting from a full geometry optimization in the ground- 
state. Rather than a state-specific geometry optimization, this allows to clearly disentangle 
geometric and electronic effects and thereby to directly compare different methods as all 
are evaluated for the same geometry (also in Fig. [s] below). The PESs shown in Fig. [2] 
are for the GGA-PBE xc functional, and we obtain essentially the same topologies for the 
three surfaces with the LDA or B3LYP. The only difference are more or less constant offsets 
between the three surfaces depending on the level of xc treatment, which is why we restrict 
the presentation for the moment to the GGA-PBE case and return to the xc discussion in 
the next section when focusing on the vertical excitation energies. 

Qualitatively, the overall obtained topology of SO, SI and S2 is perfectly consistent with 
the prevalent understanding of the azobenzene photochemistry as summarized in the pre- 
ceding section. The ground-state PES is dominated by the two metastable states, E-Ab 
and Z-Ab, separated by sizable barriers along both the rotation and inversion pathway. 
In contrast, the SI PES does not exhibit a barrier along the rotational pathway, which 
after photoexcitation of either E-Ab or Z-Ab should thus quickly lead the system to the 
well-known CI region around mid-rotation. Moreover, the Sl-FC region at E-Ab is rather 
flat, while the Sl-FC region at Z-Ab is very steep. This is perfectly consistent with the 
experimentally reported longer excited state lifetime of E-Ab compared to Z-Ab, and is also 
in line with the reported lower SI quantum yields for E— i-Z than for Z— i-E isomerization. 
Finally, the closeness of the S2 minima to the respective S0-S2 FC structures, as well as 
the separation of these minima by large barriers along both inversion and rotation suggests 
that isomerization after S2 excitation does indeed not occur on the S2 surface, but rather 
via deexcitation along CIs in other degrees of freedom than those scanned here, followed by 
motion on the SI surface. 

Essentially the same SO, SI and S2 topologies are also obtained at the RI-CC2 and 
TD-DFT level of theory, with the TD-DFT PESs obtained at different xc functional levels 
again primarily vertically shifted against each other. We present the corresponding two- 
dimensional PESs in the supplementary materiaP^and focus here instead on one-dimensional 
PES scans which allow for a more quantitative comparison. Figure |3] compiles these scans 
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CNNC dihedral angle m (deg) CNN angle a (deg) 



FIG. 3. PES scans along the rotational (left) and inversion (right) pathway, cf. Fig. [TJ Shown 
are the groundstate (SO, black), first (SI, red) and second (S2, blue) excited states, calculated each 
time with ASCF(GGA-PBE) (solid hne), TDDFT(GGA-PBE) (dotted hne) and RI-CC2 (dashed 



along the two prevalently discussed isomerization pathways, the rotational one following 
motion around the dihedral angle u and the inversion one following motion along one of 
the two CNN angles a, cf. Fig. [T} The basis of these one-dimensional PES scans are again 
optimized ground-state geometries, in which the corresponding angle was frozen and all 
other degrees of freedom were fully relaxed. 

The topological similarity, i.e. relative energetics within each PES, for the three meth- 
ods is rather striking. At the ground-state it is reflected by an almost quantitative agree- 
ment of the inversion barrier computed with DFT(GGA-PBE) and RI-CC2, 1.52 eV and 
1.74 eV, respectively. These values match also nicely with experiment and previously re- 
ported values of 1.5 eV (DFT(LDA)Pll) and 1.74eV (DFT(B3LYP)I2S1). Only slightly lower 
quantitative consistency is achieved for the SO rotation barrier, for which we compute val- 
ues of 1.81 eV and 2.05 eV at DFT(GGA-PBE) and RI-CC2 level of theory, while values 
of 1.65 eV (CASSCF/PT2P) and 2.18eV (DFT(B3LYP)IZS1) can be found in literature. At 
the RI-CC2 reference, the rotation pathway is dominated by the well known CI seam be- 
tween SO and SI state around midway rotation!- ^^ * ^^ * ^^ * ^^ *. The residual gap obtained in the 
present calculations, cf. Fig. [3| is hereby due to the fact that the geometries along the 
scan were optimized at the DFT level. For the inversion pathway RI-CC2 reveals no in- 
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tersection between SO and SI along the here displayed PES cut. Also this is in agreement 
with preceding work, which either did not find any near-degeneracies between states on the 
inversion pathway or found them only at very high energies compared to the CI seam on the 
rotation pathway In both cases, i.e. rotation and inversion pathway, ASCF(GGA-PBE) 
and TD-DFT(GGA-PBE) correctly reproduce the existence viz. non-existence of SI and 
SO state degeneracies, cf. Fig. |3j At the S2 surface, the large barrier along the inversion 
pathway is again rather well reproduced by the three methodologies, 2.08 eV (ASCF(GGA- 
PBE)), 1.55 eV (TD-DFT(GGA-PBE)), 1.50 eV (RI-CG2) when measured from E-Ab, or 
0.84eV (ASCF(GGA-PBE)), 0.82eV (TD-DFT(GGA-PBE)), 0.67 eV (RI-CC2) when mea- 
sured from Z-Ab. The same holds for the barrier along the rotation pathway on S2, where 
the corresponding values are 1.35 eV (ASCF(GGA-PBE)), 0.71 eV (TD-DFT(GGA-PBE)), 
1.03 eV (RI-CC2) when measured from E-Ab. 

We find the excellent agreement of the three methods (ASCF, TD-DFT, RI-CC2) with 
respect to the topology to also extend to other parts of the PES not contained in the hitherto 
presented scans. This is nicely demonstrated by Table [T| which compiles selected structural 
parameters of the ground-state E-Ab and Z-Ab states, as well as of those minimum energy 
structures on the SI and S2 states that are obtained after optimization from the E-Ab and Z- 
Ab FC structure. At the ground-state level we reproduce the known excellent performance 
of the DFT GGA-PBE functional in describing both geometry and relative energetics of 
both E-Ab and Z-Ab isomers as compared to both higher-level theory and experiment. 
Optimization in the SI state with FC E-Ab as starting structure yields a minimum at very 
similar geometries within ASCF(GGA-PBE), TD-DFT(GGA-PBE) and RI-CC2, which in 
turn compare nicely to the CASSCF optimized geometry reported by Cembran et al^. The 
(iNN bond length at this minimum energy geometry is essentially unchanged with respect to 
ground-state E-Ab, as one would intuitively expect for a transition depopulating the non- 
bonding n HOMO orbital. In contrast, in terms of energetics the different methodologies 
yield again a large scatter and we will return to this point in the next section. In the S2 state 
an equivalent situation is obtained. Excellent agreement in the geometries is contrasted by 
strong discrepancies in the energetics for both minimum energy structures obtained after 
optimization from the FC E-Ab and FC Z-Ab. With respect to the angular degrees of 
freedom both these minimum energy structures exhibit values in close correspondance to 
their SO counterparts. They differ largely in their strongly activated ciNN bond though. 
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TABLE I. Optimized geometry parameters of E and Z azobenzene in ground (SO) and excited 
(S1,S2) states, cf. Fig. [l]for the definition of the azo-bridge bond length c^nn and the two angles oj 
and a. Additionally shown are the relative energies IS.E of the corresponding states with respect 
to the ground-state E-Ab zero reference. None of the methods identified a stable minimum after 
optimization from SI Z-Ab, which is why the corresponding entries have been left blank in the 
table. 
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(A) (eV) 


SO 


DFT(GGA-PBE) 


180 


115 


1.26 





12 


124 


1.25 


0.58 




RI-CC2 


180 


114 


1.27 





7 


121 


1.27 


0.47 




CASSCF/PT271 


180 


115 


1.24 





4 


123 


1.24 


0.52 




ExiJ5Zl2Ql 


180 


114 


1.25 








122 


1.25 


o.dsn 


SI 


ASCF(GGA-PBE) 


180 


130 


1.25 


1.67 












TD-DFT(GGA-PBE) 


180 


131 


1.24 


1.53 












RI-CC2 


180 


128 


1.26 


2.26 












CASSCF/PT271 


180 


129 


1.25 


1.95 










S2 


ASCF(GGA-PBE) 


180 


113 


1.36 


2.53 


18 


127 


1.31 


3.60 




TD-DFT(GGA-PBE) 


180 


111 


1.34 


3.15 


30 


122 


1.31 


3.65 




RI-CC2 


180 110 /113 1.37 


4.06 












CASSCF/PT278 


180 


113 


1.35 


4.05 


8 


129 


1.29 


5.55 



This is again congruent with the expectations for a transition depopulating the 7r-bonding 
HOMO-1 orbital and supports the already mentioned interpretation that deexcitation from 
S2 occurs via CIs in other degrees of freedom than those relevant for the isomerization 
process. 

C. Vertical excitation energies 

As already mentioned in the preceding section, equivalent SO, SI and S2 topologies to 
those just discussed are obtained when using different xc functionals in the DFT, ASCF 
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FIG. 4. PES scan along the inversion pathway, calculated with ASCF (left side) and TDDFT 
(right side). Shown are the groundstate (SO, black), first (SI, red) and second (S2, blue) excited 
states, calculated each time with an LDA (dashed line), GGA-PBE (straight line) and B3LYP 
(dotted line) functional. Also shown for the TDDFT case is the pathway calculated with the 
CAM-B3LYP functional (dashed-dotted hue). 

and TD-DFT calculations. Cum grano salis corresponding PESs essentially exhibit global 
vertical shifts with respect to each other as further illustrated in Fig. |4] for the inversion 
pathway. This allows us to focus the discussion on the different levels of theory on single 
prominent points on the PES, suitably the vertical excitation energies as there experimen- 
tal data is also available as reference. Table [TT] summarizes the corresponding data. The 
benchmark against the experimental vertical excitation energies at E-Ab and Z-Ab empha- 
sizes the known high accuracy achieved by the RI-CC2 approach for the two low-lying Ab 
singlet excitationd^ and justifies its use as a theoretical reference method in our study. In 
contrast, TD-DFT and ASCF based on present-day local and semi-local functionals yield 
excitation energies that are dramatically too low, as had already been noticed in preceding 
work for this molecul d^^ * ^^ l This concerns predominantly the S2 excitation which at Z-Ab is 
underestimated by more than 1 eV. 

While thus unsatisfactory on the absolute scale, the similarity of the results produced 
by ASCF and TD-DFT as long as they are based on the same functional is notable. The 
differences are with ~ 0.4 eV largest for S2 E-Ab, while for the other three excitations 
listed in Table [TTl the two methods match to within 0.1-0.2 eV. To one end this is due 
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TABLE II. Vertical excitation energies for SI and S2 excitation at E-Ab and Z-Ab at the different 
levels of theory and from experiment. 



method 


Trans (E) 


Cis (Z) 


SO^ 


SI S2 


SI S2 




(eV) 


(eV) 


ASCF(LDA) 


2.30 2.88 


2.10 3.56 


TD-DFT(LDA) 


2.09 3.36 


2.19 3.46 


ASCF(PBE)-not spin pol. 2.27 2.75 


2.13 3.54 


ASCF(PBE) 


2.21 2.98 


2.10 3.63 


TD-DFT(PBE) 


2.15 3.39 


2.29 3.43 


ASCF(B3LYP) 


2.41 3.33 


2.30 3.83 


TD-DFT(B3LYP) 


2.53 3.63 


2.49 4.01 


TD-DFT(CAM-B3LYP) 


2.72 3.94 


2.58 4.56 


RI-CC2 


2.84 4.07 


3.00 4.51 


CASSCF/PTiZSl 


2.53 4.23 


2.72 4.49 


Exp.l™ 


2.82 4.12 


2.92 4.68 



to the sum-rule correction we employ for ASCF, which yields a spin purified state that 
is better comparable to TD-DFT^^. The plain spin-mixed ASCF approach instead yields 
excitation energies for the two states that are typically about 0.3 eV lower than the spin- 
purified ones. This would further increase the difference to the corresponding TD-DFT 
values and, worse, also to experiment. Table [TT] furthermore indicates that simple non-spin 
polarized ASCF calculationd^Sl are also in this system an alternative, effective way of tackling 
the multi-determinantal singlet problem. The corresponding values do not differ much from 
the sum- method corrected ones, and come at a significantly lower computational cost. 

To the other end the obtained similarity of ASCF and TD-DFT results is connected to 
the pronounced single-particle character of the SI and S2 excitations, which in turn is also 
the rationalization for the high accuracy of the RI-CC2 method^. This is most obvious 
for the SI state, which is essentially described by a single excitation over the entire PES 
range scanned and which is thus most straightforwardly mimicked in ASCF. However, it 

14 



also holds to some extent for the S2 state, which is of a more collective nature in the sense 
that it exhibits multiple significant TD-DFT excitation amplitudes. As apparent from the 
afore described not excessively large deviations even this can still be relatively well described 
by single-excitation ASCF due to its desirable ability to account for orbital relaxation {vide 
infra) . 

Despite the similarity, the absolute performance of both methods when based on local 
and semi-local functionals against RI-CC2 and experiment is still a concern. Particularly the 
deviation in the S2 Z-Ab vertical excitation energy exceeds the one commonly found for low- 
lying singlet excitations in organic molecule^^^**^. Visual inspection of the involved frontier 
orbitals, cf. Fig. |6] below, suggests that this might be related to some charge-transfer (CT) 
aspect of the excitations, which in particular for S2 around Z-Ab shifts charge between the 
central azo-group and the phenyl-moieties. This interpretation receives some quantitative 
support by an evaluation of Tozer's CT A-parametei'^. Measuring the spatial overlap in a 
given excitation, A values towards unity indicate that occupied and virtual orbitals involved 
in the excitation occupy increasingly similar regions of space. In contrast to Rydberg and CT 
excitations such " short-range excitations" should then be much better amenable to TD-DFT 
based on local or semi-local xc functional ij^^ * ^^ * ^^ -*. For the present case, it is indeed precisely 
for S2 and around Z-Ab that we obtain the lowest A- values around 0.5, while for SI and for 
S2 towards E-Ab A- values lie consistently around 0.7 or higher. Further support for the CT 
picture comes then also from the larger reduction of the underestimation of S2 at Z-Ab when 
going to the hybrid functional level. Whereas for SI and S2 at E-Ab the admixture of exact 
exchange reduces the error with respect to the RI-CC2 from ~ 0.7eV to ~ 0.4 eV, at Z-Ab 
the S2 error goes from the larger ~ 1.1 eV equally down to ~ 0.5 eV, cf. Table [TTj At the level 
of the Coulomb-attenuated CAM-B3LYP functional this remaining deviation is then further 
reduced to the order of 0.1 eV throughout. Also this is consistent with the interpretation 
of some overall CT character of the low-lying singlet excitations, which this functional with 
its varying degree of exact exchange at short and long range is specifically supposed to 
tackl^^^'^^. At this functional level the agreement reached with respect to experiment is 
thus essentially en par with that of the correlated wave function reference techniques, and 
we speculate that the same would approximately hold for CAM-B3LYP-based ASCF (which 
unfortunately is presently not available to us). 
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FIG. 5. Scan along the inversion pathway, comparing the S2 PES as calculated with TD- 
DFT(GGA-PBE) when following the HOMO-l-^LUMO excitation (red, dotted line) and when 
following the transition with predominant vr — )• vr* character (green dotted line). Only the latter 
approach yields the correct PES topology with sizable barrier around mid-inversion as compared 
to RI-CC2 (dashed line) and ASCF(GGA-PBE) (sohd line). Also shown are the CT A-valuePl 
for the two TD-DFT curves (see text), as well as the LUMO and the "wrong" HOMO-1 orbital at 
mid-inversion. The latter nicely reveal the pronounced CT-character of this transition. 

D. TD-DFT and inversion path barrier 



The good agreement of ASCF, TD-DFT and RI-CC2 in terms of overall PES topology 
particularly around mid-inversion might come as a bit of a surprise in view of earlier studies 



molecule^^^^J*^. In such cases deviations between ASCF and TD-DFT are often readily 
attributed to the "simplicity" of the prior theory. Alternatively, "collective" character of 
an excitation as judged from the existence of several significant TD-DFT amplitudes is 
also cited as reason for the failure of "single excitation" restricted ASCF. In turn, when it 
comes to differences between TD-DFT and higher level wave function theories, the deficiency 
of semi-local TD-DFT to deal with CT-character of excitations is a frequently encountered 



that reported significant discrepancies for thid^ or for vr-bond twisting paths of comparable 
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FIG. 6. Energetic positions of GGA-PBE KS frontier orbitals along the inversion pathway as 
resulting from a self-consistent ground-state calculation (upper panel) and as resulting from a 
self-consistent ASCF calculation for the S2 excitation (lower panel). Additionally shown are the 
corresponding KS ground-state orbital shapes at Z-Ab (left) and E-Ab (right). 

rationalization. Instead, in the present case yet another difficulty of TD-DFT applies, namely 
state-crossings, and we find it instructive to point out that this can easily lead to wrong 
assessments, in particular as we find the allegedly "simpler" theory ASCF to be significantly 
more robust with respect to this issue. 

The reason why a similar S2 topology with sizable barrier around mid-inversion was found 
for ASCF, TD-DFT and RI-CC2 in Fig. [3] is that we do not simply plot the values following 
the transition from the second highest occupied to the lowest unoccupied orbital without 
considering the orbital character in the case of TD-DFT. Instead, we always specifically 
follow the excitation that exhibits the largest amplitude for the targeted transition between 
the IT and tt* KS orbitals. As shown in Fig. |5] the two procedures do not lead to the same 
result around mid-inversion, and only the approach that tracks the correct transition yields 
a PES topology that is in agreement with both ASCF and the reference RI-CC2 data. In 
contrast, the approach that merely monitors the HOMO-1 to LUMO TD-DFT excitation 
yields a very wide plateau-region along the inversion path, precisely as reported previously 
by Tiago et al^ for this system. 

The source for this difference is clearly apparent from Fig. |6} which shows the evolution 
of the energetic position of the GGA-PBE KS frontier orbitals along the inversion pathway. 
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At both E-Ab and Z-Ab these orbitals exhibit the ordering HOMO-1 (vr), HOMO (n) and 
LUMO (tt*) as intuitively expected from the n ^ n* (HOMO to LUMO) and n ^ ir* 
(HOMO-1 to LUMO) character of SI and S2, respectively. However, around mid-inversion 
the ground-state KS orbital energies of the states that at E-Ab and Z-Ab correspond to 
HOMO-1 and HOMO-2 cross. Inspection of the TD-DFT amplitudes reveals that the correct 
TT — 7- TT* TD-DFT excitation then corresponds to a transition primarily between this HOMO- 
2 level and the LUMO. 

Here, it is intriguing to note that also in Hartree-Fock as the basis for the RI-CC2 
reference, the frontier orbital ordering does not correspond to the aforementioned vr, n, and 
TT* sequence. Still, the two lowest energy RI-CC2 excitations anywhere on the PES parts 
scanned in this study have predominant amplitude just exactly for the transitions expected, 
i.e. n — i- TT* for SI and tt — t- vr* for S2. As such we interpret the TD-DFT result which 
places the vr — ?■ vr* transition at a higher lying excitation than the second transition around 
mid-inversion as the inability of linear-response theory to cope with the bad ground-state 
orbitals offered by GGA-PBE. This is then also consistent with the observation that in TD- 
DFT the S2 excitation typically exhibits larger amplitudes for more than one single-particle 
transitioii^. 

Intriguingly, the allegedly simpler theory ASCF is much less affected by this limitation as 
it allows for orbital relaxation under the excitation constraint. Figure |6] demonstrates that 
the self-consistent orbitals obtained under the S2 population constraint no longer exhibit any 
state crossing along the inversion pathway. Once self-consistency is achieved, the constraint 
of a depopulated HOMO-1 and a populated LUMO leads always to a situation where the 
HOMO-1 corresponds to the tt state and the LUMO to the vr* state, i.e. the computed 
excitation energy corresponds indeed exactly to the transition we wanted to model. We found 
this type of robustness to hold for both SI and S2, everywhere on the PES parts we scanned, 
and for whatever xc functional we used. Especially for dynamical simulations or mappings 
of larger parts of the PESs this is in our eyes an important asset. It is also particularly 
remarkable as at the hybrid functional level (B3LYP and CAM-B3LYP) the ordering of the 
ground-state KS levels differs from expected vr, n, vr* sequence essentially everywhere on the 
PES parts we scanned. This made it rather cumbersome to track the correct transitions in 
TD-DFT, which, however, was the prerequisite to obtain the consistent agreement of the TD- 
DFT PES topologies with respect to ASCF and RI-CC2 reported above. For corresponding 
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method comparisons the ease with which unidentified state-crossings can impair the TD- 
DFT results is hereby particularly consequential, as it may readily lead to wrong assessments. 
Discrepancies between TD-DFT and ASCF that in reality are due to an unidentified state- 
crossing in TD-DFT may lead to the dismissal of the allegedly "simpler" ASCF theory^. As 
shown in Fig. [5] the switch of the excitation character induced by the state-crossing around 
mid-inversion gives furthermore rise to small A values for the wrong TD-DFT transition. 
This bears the danger to assign the discrepancy in the PES topology of this transition with 
respect to the RI-CC2 reference incorrectly to the deficiency of present-day functionals in 
describing CT excitations. 

IV. CONCLUSION AND OUTLOOK TO SURFACE-MOUNTED 
AZOBENZENE 

In summary we have systematically computed ground- and low-lying singlet excited state 
PESs that are of relevance for the isomerization dynamics of gas-phase azobenzene. Our 
results demonstrate that sum-method corrected ASCF yields global PES topologies, i.e. 
relative energetics within one PES, that agree very well with those of TD-DFT at the 
same xc functional level and with accurate RI-CC2 reference data. Previous contradictory 
reports concerning the agreement of ASCF and TD-DFT suffered from unresolved state 
crossings in the TD-DFT calculation^, while the orbital relaxation possible in ASCF makes 
this approach very robust with respect to this issue. The now unanimously obtained PES 
topologies of SO, SI {n — )■ vr*) and S2 (vr tt*) states are furthermore quite consistent with 
existing experimental data concerning the photo-isomerization mechanism. 

When based on the same xc functional sum-method corrected ASCF and TD-DFT agree 
to within 0.1-0.2 eV for the SI state that is most relevant for the isomerization of free 
azobenzene, and to within 0.4 eV for the S2 state. This suggests ASCF as a promising 
route to larger azobenzene-containing systems, where TD-DFT becomes computationally 
untractable. This concerns predominantly surface-mounted azobenzene, where preliminary 
calculations for the adsorption at coinage metals indicate that the low lying excited states 
largely retain their molecular charactei^. Particularly appealing in this context is that we 
obtain the correct PES topology for gas-phase azobenzene already at the level of semi-local 
xc functionals, which are presently the unbeaten workhorse for metal adsorption studies. 
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Notwithstanding, at this level of xc treatment the vertical excitation energies produced 
by ASCF and TD-DFT are grossly underestimated. Our analysis attributes this primarily to 
some charge-transfer character of the SI and S2 excitations, which the local functionals are 
unable to grasp. As to be expected hybrid (B3LYP) and even more so long-range corrected 
hybrid (CAM-B3LYP) functionals improve on this situation. Conserving the overall topol- 
ogy they primarily induce global upward shifts of the excited state PESs, leading to vertical 
excitation energies at the CAM-B3LYP level that are roughly en par to the correlated wave 
function reference techniques. 

While the reduction of self-interaction error achieved with the advanced functionals is 
therefore of paramount importance to adequately describe the molecular excitations, one has 
to recognize that mere admixture of exact exchange does not seem to be the right pathway 
for adsorption at metal surfaces (at least when judged from the few seminal ground-state 
studies performed to date^^^^^^. There, semi-local functionals are still more or less the only 
tractable approach. For future applications to surface-mounted azobenzene this dictates a 
cautious approach carefully assessing what can and what cannot be addressed with semi- 
local functional based ASCF. The present results give good hope that the qualitative PES 
topology may be retrieved, with recent methodological ASCF extension^^Mo] particular 
interest in case of appreciable adsorbate-surface hybridization. In the best of cases this is 
already sufficient to conclude on the isomerization mechanism. As to the absolute excitation 
energies the situation is less clear. Semi-local levels of xc treatment that are adequate for the 
extended metallic states will be insufficient for the localized molecular states. We speculate 
that an insightful mixture of extended supercell surface calculations and suitable molecular 
model system reference calculations might nevertheless yield highly valuable insight to this 
end. 
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